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ABSTRACT 


A mathematical method from random variable theory was used to 
analyze simultaneous variation of engine performance and stage inert weight 
for some typical 2-stage rocket vehicles. This represents an extension of 
the conventional sensitivity analysis which is often performed for conceptual 
designs. The study was incomplete in that too few examples were analyzed 
for conclusions to be drawn regarding the usefulness of the approach, but 
mathematical consistency was shown. 
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DEFINITION OF SYMBOLS 


NOMENCLATURE 


g Gravitational acceleration (9. 805 m/sec 2 ) 

I Propulsive specific impulse, seconds 

M Mass, kg 

p (x) Density function for x 

P (x) Distribution function for x 

v Velocity, m/sec 

y Dummy variable 

z Dummy variable 

X Stage propellant fraction 

ju . Vehicle mass ratio, (mass fully loaded)/ (mass with impulse 

propellant expended) 

£ Mass ratio, Mt/Ml 

SUBSCRIPTS 

L Payload 

T . , Total 

1 First stage 

2 Second stage 


Where appended to functional relationships (e. g. , g lt g 2 , etc. ) , numerical 
subscripts do not refer to stages. 
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PARAMETRIC SENSITIVITY ANALYSIS 
FOR STAGED PROPULSIVE VEHICLES 

SUMMARY 


A mathematical method from random variable theory was used to analyze 
simultaneous variation of engine performance and stage inert weight for some 
typical 2-stage rocket vehicles. This represents an extension of the conven- 
tional sensitivity analysis which is often performed for conceptual designs. 

The study was incomplete in that too few examples were analyzed for conclu- 
sions to be drawn regarding the usefulness of the approach, but mathematical 
consistency was shown. 

BACKGROUND AND INTRODUCTION 


During preliminary and conceptual design phases of studies of flight 

systems, it is common practice to perform sensitivity analyses to determine 

the effect of variation of system parameters that cannot be precisely predicted 

in the course of such early studies. These sensitivity data are commonly in 

the form of derivatives such as change in vehicle gross weight as a function of 

change in propulsion I of one of the stages, and change in system gross 

sp 

weight versus the sub-parameters of interest. In a plotted form one can 
readily see, for example, how much the system gross weight will be affected 
by a 1% variation in second stage inert weight. It is of interest to investigate 
the combined effect of simultaneous variation in more than one parameter. 

For example: What is the combined effect of uncertainty in delivered I and 
inert weights in a i- or 2-stage vehicle? sp 

A method that can be used to simulate simultaneous variation in critical 
sub-parameters is the treatment of these sub-parameters as random variables. 
This method can be readily employed whenever it is possible to express the 
desired vehicle performance parameters in mathematical form. Gross weight- 
to-payload ratio may be expressed as a mathematical function of the various 
sub-parameters such that the pertinent partial derivatives may be analytically 
evaluated. This may be done with a simple ideal velocity model of rocket 
vehicle performance. Even this simple model will provide valuable insight 
into the likelihood of achieving design goals and into the selection of certain 
parameters, such as staging velocity, to minimize the development risk. 



The study reported herein is incomplete in that too few examples were 
analyzed to achieve a good evaluation of the usefulness of the technique. How- 
ever, a mathematical approach was formulated using the random variable 
theory, and enough calculations were made on a digital computer to show that 
the approach is at least mathematically consistent. 


PURPOSE 


The purpose of this internal note is to report on an exploratory in- 
vestigation of a method of analysis for treatment of vehicle performance by 
random variable techniques and to give example results of the analysis for 
representative rocket propelled vehicles. 

Application of the random variable theory as suggested above requires 
some information about the variables considered to be random. This infor- 
mation is normally in the form of a distribution density function. Such a 
function must be either known or assumed. 


As a simple example, suppose that it is desired to determine the density 
function for required mass ratio (initial weight/cutoff weight) , knowing the 
density function for engine performance in terms of I . The density function 

of a random variable is defined as the derivative of the distribution function 
of the variable, which is the probability of the variable being smaller than a 
given value. This relation is sketched graphically below, where 



a> 
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Where an assumed function p [j] is provided, with an arbitrary ordinate, 
it must be normalized by the relation that follows from the definition of the 
distribution function P [l] and equation (1): 

CO 

f 0 P [II dl = P[») =i (2) 

There is a distribution function for the required mass ratio, /x, that 
corresponds to the distribution function for I. For a given value I' there is 
a corresponding value fx\ The probability that (x is greater than ju' is equal 
to the probability that I is smaller than I f , i, e. , 


P [m* ] " 1-P [I 1 ] 

Usually it is the density function that is desired. 


(3) 


P t M ] 


bp Pm 1 

0 IX 


Substituting from (3), 

P [Ml = 


-op[i ] 

d(X 

-3P[l] 


81 


= PH 


Ml 

81 


81 

8jU 


The relation between fx and I is 

AV 


H = exp 


“gl“ 


where AV is the ideal AV (disregarding gravity and drag losses) . 


dix _ -AV AV 

81 gl 2 GXp gl 


(4) 


(5) 


( 6 ) 


(7) 
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Substituting (7) and (5) allows determination of p [pi]. 

Analysis of a vehicle could proceed using X as a random variable 
instead of pi. However, it may be interesting to know the density function 
of pt. It will be assumed in the following analysis that a preliminary step 
has determined p [pi ] . 

The equation for total mass of the second stage of a 13-stage vehicle or 
the total mass of a single stage vehicle is 



The joint density function [ pi , A] is known: 


PpiA ] = P ] P [A] (11) 

This follows from the assumption that p [pi] and p [A] are independent. 

For convenience, letting £ - Mt/Ml, the joint density function . [ A] 
if then found by: f 

_9£i Mu 

Py, r (y>£) “ P/z,,\[gi (y.£) ( §k (y»£) J M Qm (12) 

Ml Mi 

8 £ 9 £ 
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g ( izxi 

px£ (y,g) = p ix fa) px\ y <*-£) 


JJL 

ax 

8 ft 

dn 


ax 

9£ 

9£ 


( 12 a) 


8 X 

0^t 

0 /i 

ax 

0£ 


and 9 y are both zero. 

dt 


= 1 

JL 

sg 


£ . U ry). 


y ( i- a 


( 12 a) becomes 

p[y>£3 = [M] EX 

It is assumed that p [y,£] ~ p [y] p [£]. 


i-y 


y (i-t)' 


6 (i-y) 


L* (i-t)J 



,o 

l-y 


o y (l-g) 2 


Therefore, p [£] may be found by integrating over y: 

CO 

P (£) = f P(y,£) dy 

*_nn 


(13) 


(14) 


(15) 


Since y - (i this is equivalent to 

(i max. 

P(S > = .Uin. P y ’ £ < p ' *1 dp 


(15a) 


This integral may be performed for a number of values in the range of 
£. If it is also desired to get the distribution function for £, this may be found 

b y 


p (£) = f P (£’) d£j 

£min. 


(16) 


For a 2-stage vehicle, £ in equation (16; becomes £ 2 311(1 £ - £ 2 , 

since the second stage is the payload of the first stage. £ 2 is now a random 
variable defined by equation (16). Analagous to equation (8) , 
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y _ £? ^-9 M 9 

l-(l**l)Ml 

Two dummy variables are required: 


( 17 ) 


yi = Mi 


Z 1 = ?2 


The same type of inverse relations hold: 

Mi = gi (yi, z if £ ) 3 y t 

f2“ (yi. zi, £ ) = zi 

*■1 = g3 (yi. z ls & ) 3 


JJlzZi) 
yi U-n 


The joint density function is now: 

P [M, *i. k) = P [Mi] P M p [kl 
and the joint density function p r [y 1# z l9 £ ] 

y 1 z * £ 


is then founded by: 


P. 


y,z, £ 


[yi» z i> £ ] = PmH 2 [gi, g 2 . gs] 


All of the partial derivatives are zero except 

_9&l 


_9gi - _9g2 = j_ and 


i-yi 


(18) 


(19) 


( 20 ) 


Jgl 

9go 

9ga 

3M 1 

3Mi 

9Mi 

Jgl 

3 go 

Pga 

3z t 

3zj 

8zj 

9gi 

-3&L 

_Sga_ 

3£ 

3£ 

9£ 


( 21 ) 


8M1 9 z i ~ 3£ 

p [£] may be found by integrating yj and Zj 


yi (1- n 


( 22 ) 


ru n _ fk max. rn< max. . . .. . 

P[fl - Vmin. i/mta. P (l ' £ * £ > dp ' 


(23) 
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The distribution function P [£] is again found by equation (16). 


RESULTS 

Explanation of Terminology 


For the non-technical reader, an explanation is in order regarding 
the format of data to be shown. First, what is a density function? The 
density function represents something that is not precisely known and about 
which there is some uncertainty; it does, however, provide information on 
probability. If we have a function whose exact value is not known, and we know 
that there is a 25% probability that its value is between 3 and 4, we can draw 
a histogram, as shown below. 



If we have more information about this function, that is, if we know 
the probability of its being between 1 and If and 2, etc. , we can draw a 
histogram showing more detail, as in the following sketch. 
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If we divide the region into arbitrarily small elements instead of a 
stepped curve, then a continuous curve results. This is a density function. 

If the density function is normalized, any point on the curve represents the 
probability that the value of the function will fall within a unit range about the 
point. For example, consider a density function plotted versus a specific 
impulse obtainable from a rocket engine; a point plotted at a value such as 400 
seconds represents the probability that the I S p obtained will lie between 399. 5 

and 400. 5 seconds. The importance of normalization lies in the fact that the 
integral from zero to infinity of a normalized density function is exactly equal 
to 1. 


This integral then forms what is known as the distribution function. The 
distribution function is either constant or monotonically increasing and re- 
presents the probability that a random variable will be less than the abscissa 
value against which the distribution is plotted. If at an I g p value of 400 the 

distribution function is equal to 0. 4, then the probability is 40% that the I S p 

will be less than 400, and 60% that it will be greater than 400. 

It is convenient for the analyst using this random variable technique 
to make an "eyeball" estimate of the uncertainties in his sub -parameters in 
the form of density function. These may then be employed to show the effects 
of variations in sub-parameters on overall system performance parameters. 

It is convenient that the analyst should not have to produce a normalized density 
function. An ordinate scale of 0 to 5 is typically convenient. Normalization 
can be readily carried out by the computer, prior to its using the input density 
functions in the calculations. 


Example Calculations 

Computation of p [£ ] and P [f ] was performed by a digital computer 
routine for certain 2-stage vehicle examples. Inputs were assumed density 
functions for first and second stage I g p and propellant fractions, payload mass, 

and overall impulsive velocity increment for the vehicle, and staging velocity. 
Input density functions were given an arbitrary scale of 0 to 5 and were normali- 
zed by the computer routine. The following examples were given: a lunar 
direct flight and return vehicle, a Saturn-type 2-stage to orbit expendable 
vehicle, and a reusable space transport vehicle with rocket propulsion. Results 
from the calculations were automatically plotted on the SC-4020 plotter. 


8 



For the lunar direct flight vehicle a single-stage lander and single- 
stage return vehicle staging necessarily occurs at the lunar surface for the 
case shown. The given impulsive velocity increments are not typical of those 
currently used for lunar mission analysis, but they serve to illustrate typical 
results of tills analysis. Both stages were assumed to have the same engine 
performance. The landing stage is assumed to have abetter propellant fraction 
because it is larger and because it is not required to be well insulated for 
extended storage of cryogenic propellants. Figures 1 through 5 give the 
results for this vehicle. Irregularities in Figure 2 result from the finite 
difference nature of the numerical integrations employed. 

The next two examples, where the staging velocity is not fixed by mission 
characteristics, were calculated for several staging velocities. This allows 
selection of a staging velocity that will give reasonable all-up weights with 
minimum sensitivity to uncertainty in inert weights and engine performance. 
Accumulated errors in numerical integrations may be noted by the fact that 
distribution function curves for these cases do not terminate exactly at i. 0 as 
they should. 

The random variable approach is not limited to the relatively simple 
cases shown. It could be extended to more than two stages, and variation 
in nominal propellant fraction with propellant loading could be incorporated 
without much difficulty. 



INPUT DENSITY FUNCTION INPUT DENSITY FUNCTION 


Lunar Direct Flight Vehicle 


One-stage Cryogenic Lander 
Return Payload Mass - 5670 Kg 
Return Ideal Delta V = 2. 75Km/Sec 
Landing Ideal Delta V 3 3. 697 Kxn/Sec 
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FIGURE 1 
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DENSITY FUNCTION FOR MASS 


Lunar Direct Flight Vehicle 
One-stage Cryogenic Lander 
Return Payload Mass = 5670 Kg 
Return Ideal Delta V - 2. 75Km/Sec 
Landing Ideal Delta V - 3, 697 Km/Sec 



FIGURE 2 
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DISTRIBUTION FUNCTION FOR" MASS 


Lunar Direct Flight Vehicle 

One -Stage Cryogenic Lander 
Return Payload Mass = 5670 Kg 
Return Ideal Delta V - 2. 75Km/Sec 
Landing Ideal Delta V * 3. 697 Km/Sec 



Return Stage Initial Mass, Kg 


FIGURE 3 



DENSITY FUNCTION FOR MASS 


Lunar Direct Flight Vehicle 

One-Stage Cryogenic Lander 
Return Payload Mass = 5670 Kg 
Return Ideal Delta V = 2. 75Km/Sec 
Landing Ideal Delta V = 3. 697 Km/Sec 



FIGURE 4 
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DISTRIBUTION FUNCTION FOR MASS 


Lunar Direct Flight Vehicle 
One-Stage Cryogenic Lander 
Return Payload Mass - 5670 Kg 
Return Ideal Delta V- 2. 75Km/Sec 
Landing Ideal Delta V- 3. 697 Km/ Sec 



FIGURE 5 
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INPUT DENSITY FUNCTION INPUT DENSITY FUNCTION 


Saturn- Typo Two Stage to Orbit Launch Vehicle 

Payload Mass = 125 Metric Tons 

Total Ideal Velocity Increment - 9150 M/Sec 
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Specific Impulse, LB-Sec/LB 
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•Stage Propellant Fraction 
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Specific Impulse, LB-Sec/LB 



FIGURE 6 
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DENSITY FUNCTION FOR MASS 


Saturn- Type Two-Stage-to-Orbit Launch Vehicle 

Payload Mass * 125 Metric Tons 

Total Ideal Velocity Increment - 9150 M/Sec 



6.000X10+ 05 8. OOOXlO+05 1.000X10' 

Second Stage Initial Mass, Kg 


FIGURE 7 


Saturn- Type Two Stage to Orbit Launch Vehicle 

Payload Mass - 125 Metric Tons 

Total Ideal Velocity Increment = 9150 M/Sec 



6. OOOXIO* 05 8. OOOX10 + ® 5 1.000X10 + ° 5 

Second Stage Initial Mass, Kg. 

FIGURE 8 
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DENSITY FUNCTION _ JR MASS 


Saturn-Type Two Stage to Orbit Launch Vehicle 

Payload Maas * 125 Metric Tons 

Total Ideal Velocity Increment * 9150 M/Sec 



FIGURE 9 
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DISTRIBUTION FUNCTION FOB MASS 


Saturn- 'type Two Stage to Orbit Launch Vehicle 

Payload Mass « 125 Metric Tons 

Total Ideal Velocity Increment = 9150 M/Sec 



FIGURE 10 
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INPUT DENSITY FUNCTION INPUT DENSITY FUNCTION 


Reusable Space Transport, 2 Stage 
F-l First Stage, Hg-3 Second Stage 
Payload Mass - 11.35 Metric Tons 
Total Ideal Velocity Increment = 8. 85 Km/ Sec 



Specific Impulss, LB-Sec/LB 



Stage Propellant Fraction 


First Stage l 
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1 
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284 285 286 287 288 289 290 291 
Specific Impulse, LB-Sec/LB 



Stage Propellant Fraction 


FIGURE 11 
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Reusable Space Transport, 2 Stage 

F-l First Stage, Hg-3 Second Stage 

Payload Mass = 11. 33 Metric Tons 

Total Ideal Velocity Increment = 8, 85 Km/Sec 


4. 000X10 •° B 



5. 000X10+® 1. 000X10+* 

Second Stage Initial Mass , Kg 


FIGURE 12 



DISTRIBUTION FUNCTION FOR MASS 


Reusable Space Transport, 2 Stage 

F-l First Stage, Hg-3 Second Stage 
PayLoad Mass = 11,35 Metric Tons 
Total Ideal Velocity Increment ~ 8. 85 Km/ Sec 

1.2 



FIGURE 13 
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DENSITY FUNCTION FOR MASS 


Reusable Space Transport, 2 Stage 

F-i First Stage, Hg~3 Second Stage 

Payload Mass “ ii. 35 Metric Tons 

Total Ideal Velocity Increment = 8. 85 Km/ Sec 



FIGURE 14 
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DISTRIBUTION FUNCTION FOR MASS 


Reusable Space Transport, 2 Stage 

F-l First Stage, Hg-3 Second Stage 

Payload Mass = 11. 35 Metric Tons 

Total Idee! Velocity Increment = 8. 85 Km/ Sec 
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